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Abstract. The ZZ Ceti star KUV 02464+3239 was observed over a whole season at the mountain station of Konkoly 
Observatory. A rigorous frequency analysis revealed 6 certain periods between 619 and 1250 seconds, with no shorter period 
modes present. We use the observed periods, published effective temperature and surface gravity, along with the model grid 
code of Bischoff-Kim, Montgomery and Winget 1 1] to perform a seismological analysis. We find acceptable model fits with 
masses between 0.60 and 0.70 Mr.,. The hydrogen layer mass of the acceptable models are almost always between 10 -4 and 
10~ 6 M t . In addition to our seismological results, we also show our analysis of individual light curve segments. Considering 
the non-sinusoidal shape of the light curve and the Fourier spectra of segments showing large amplitude variations, the 
importance of non-linear effects in the pulsation is clearly seen. 
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INTRODUCTION 

The luminosity variations of KUV 02464+3239 have 
been discovered by Fontaine et al. [6]. They identified 
it as a DAV pulsator using an about 50-minute-long ob- 
servation made in 1999. In this paper we present some 
results of the Fourier and seismological analyses of the 
star based on our observations. We also present our in- 
vestigations of selected light curve segments. 

Based on its atmospheric parameters (77 e ff = 11 290 K, 
logg = 8.08, l6J) the star is located near the red edge 
of the DAV instability strip. Non-sinusoidal light curve 
shapes, short-term amplitude variations and the presence 
of linear combination peaks in the Fourier spectra char- 
acterize the pulsation of similar cool DAV stars. Investi- 
gations of these complex pulsators can allow of to know 
more about the dynamics behind the observed light vari- 
ations. 



OBSERVATIONS 

We started white light observations of the star in Octo- 
ber 2006 at Piszkesteto (the mountain station of Konkoly 
Observatory) with the lm RCC telescope and a Prince- 
ton Instruments VersArray:1300B CCD camera. The ob- 
servations ended February 2007. Covering the whole 
observing season we obtained data on 20 nights. Our 
longest time string obtained on JD 2454068 covers 
11.2h. 



We used standard IRAF0 routines during the photo- 
metric reduction process. Data analyses were made by 
use of the MuFrAn (Multi-Frequency Analyzer) pack- 
age HO] and the time series analysis program Pe- 
riod04 lHHl . Frequency values are derived in cycle/day 
(c/d) units. 

KUV 02464+3239 was found to be a multiperiodic 
pulsator showing large amplitude and long period pul- 
sation modes in accordance with its position in the DAV 
instability strip. Another characteristic of the light curve 
is its strongly non-sinusoidal shape. This means that non- 
linear effects play an important role in the observed vari- 
ability. 



RESULTS OF FOURIER ANALYSES 

Since we observed KUV 02464+3239 through a whole 
season, we had the opportunity to form groups of nights 
and analyse them independently. We selected four sub- 
sets of nights: four nights from October (Subset 1), the 
longest time string (Subset 2) from November and a one- 
and four-night dataset from December (Subset 3 and 4, 
respectively). Table[T]shows the data subsets' log of ob- 
servations. The Fourier transforms of the subsets can be 



1 IRAF is distributed by the National Optical Astronomy Observato- 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Sci- 
ence Foundation. 



TABLE 1. Log of observations of KUV 02464+3239 data 
subsets. N denotes the number of data points, 8T shows the 
length of data sets. 



Subset 


Date 


Start time 


N 


ST 


No. 


[UT] 


[HJD-2450000] 




[h] 


1 


2006 Oct 06 


4014.577 


210 


2.02 


1 


2006 Oct 07 


4015.584 


117 


1.31 


1 


2006 Oct 09 


4017.546 


216 


2.37 


1 


2006 Oct 1 1 


4019.543 


264 


2.98 
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2006 Nov 28 


4068.206 


1163 


11.17 


3 


2006 Dec 1 1 


4081.169 


732 


8.72 


4 


2006 Dec 14 


4084.170 


1048 


10.92 
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2006 Dec 15 


4085.179 
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10.56 
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2006 Dec 16 


4086.187 
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2006 Dec 19 


4089.208 
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6.90 
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FIGURE 1. Fourier spectra of data subsets. Numbers of sub- 
sets and corresponding window functions are inserted in the 
upper-left and right corners of each panel, respectively. We de- 
noted by arrows the 6 frequencies we accepted as normal modes 
in the last panel. 

seen in Fig.Q] 

In the panels of Fig. Q] we can follow the remark- 
able variations in the amplitudes of pulsational peaks. 
For example, the amplitudes of peaks around 100 c/d de- 
creased strongly, while the peak at ^87 c/d became dom- 
inant from Subset 2 to Subset 4. Amplitude variations 
can made the finding of pulsational modes difficult in the 
Fourier spectrum of the whole light curve. In such case 
the Fourier transform (FT) of the whole dataset gives 
an averaged amplitude spectrum. We performed detailed 
analyses of the data subsets and accepted 6 frequencies 
as normal modes of pulsation. We denoted these fre- 
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FIGURE 2. Fourier spectrum of the whole light curve. The 
6 accepted normal modes are denoted by arrows. The window 
function is given in the insert. 



TABLE 2. Frequency, period and amplitude values of 
the 6 accepted modes derived by the whole light curve. 





Frequency 


Period 


Ampl. 




[c/d] 


[MHz] 


[s] 


[mmag] 


/l 


69.1060(3) 


799.838 


1250.253 


4.4 


h 


86.9879(1) 


1006.804 


993.242 


13.2 


h 


99.7516(1) 


1154.532 


866.151 


9.5 


h 


104.2617(1) 


1206.733 


828.684 


11.6 


h 


111.1095(2) 


1285.989 


777.611 


5.5 


k 


139.5157(3) 


1614.765 


619.285 


4.0 



quency values by arrows in the last panel of Fig. Q] At 
these frequencies we found at least in one subset a very 
significant and well-determined peak. 

Fig. |2] shows the Fourier spectrum of the whole 
dataset. We marked by arrows again the frequencies 
we accepted utilizing the analyses of data subsets. Fre- 
quency, period and amplitude values derived by the 
whole light curve are listed in Table [2] All of the modes 
can be found in the long period regime, between 619 and 
1250 s. After pre-whitening with these frequencies we 
determined additional 7 in the residual spectrum. How- 
ever, they have lower probability to be normal modes. 
Most of them are close to the formerly accepted ones 
and one of them is the first harmonic of the dominant 
(J2) mode. We decided to accept the 6-frequency light 
curve solution - which is supported by the analyses of 
subsets - and not to use all the 13 (or 12 without the har- 
monic mode) frequencies for the seismological analysis 
of KUV 02464+3239. 



MODELING 

We built a grid of stellar structure models considering the 
effective temperature and surface gravity values deter- 
mined by spectroscopy [6] and their estimated uncertain- 
ties (a) [7]. Since we need mass parameters for modeling 
instead of logg values, we used seismological masses of 
DA stars determined by Bradley (J] to estimate the mass 
range we have to cover with the grid. We varied the effec- 
tive temperature and mass parameters between 10 800 - 



TABLE 3. Model paremeters of the 6 selected models and 
the corresponding r.m.s. values derived from the observed 
and model periods. Models having 1=1 solutions for the 3 
largest amplitude modes are typeset in boldface. 



M* (logg) 


T eff 


-logM H 


x 




^ r.m.s. 


[Mq] 


[K] 








[s] 


0.615 (8.03) 


11800 


4.0 


0.7 


0.3 


1.51 


0.625 (8.04) 


11000 


7.4 


0.5 


0.2 


1.50 


0.645 (8.07) 


11400 


5.2 


0.9 


0.2 


1.33 


0.650 (8.08) 


11800 


4.6 


0.6 


0.1 


0.93 


0.680 (8.13) 


11800 


5.0 


0.5 


0.1 


1.26 


0.685 (8.14) 


11400 


4.8 


0.6 


0.4 


1.12 



1 1 800 K (-1 1 290±500 K) and 0.525 - 0.74 M* (logg = 
8. 08 ±0.1 dex), respectively. The mass of the hydrogen 
layer was changed between 10~ 4 — 10 _8 M* and we fixed 
the mass of the helium layer at lO^M*. While we can 
not exclude thinner helium layers, based on a second 
scan we found that with this value we get better solu- 
tions. We varied the core parameters X Q (the central oxy- 
gen abundance) and Xf m (the fractional mass point where 
the oxygen abundance starts dropping) between 0.5 - 0.9 
and 0.1 - 0.5, respectively. Step sizes were 200 K (T e g), 
0.005 M (AT,), 10-°- 2 M* (M H ) and 0.1 (X andX fm ). 

We used the White Dwarf Evolution Code (WDEC) 
last modified by Bischoff-Kim, Montgomery and Winget 
(jl to build our model grid and get pulsation periods of 
a model. Detailed descriptions of the code can be found 
in Lamb & van Horn MlOll and Wood B15I1 . We applied 
the integrated form of the WDEC developed by Metcalfe 
lfl2tl . To find model solutions with periods closest to 
our observed ones we used the fitting routine f it per 
written by Kim [8J] . 

Our main criteria to select a model as an acceptable 
fit were a low r.m.s. value calculated from the observed 
and model periods and at least 3 / = 1 solutions for the 
6 modes. We let all 6 modes be / = 1 or 2 for the fitting 
procedure, but assuming better visibility of / = 1 modes, 
we preferred models that give at least 3 I — 1 solutions. 
Table [3] shows our 6 selected models. They fulfil these 
two criteria and have stellar masses within the 1 rj logg 
range. While the 0.625 M Q model has a thinner hydrogen 
layer, the other 5 models have hydrogen layer masses 
between 10~ 4 -6 • 10~ 6 M*. 

By further investigations of the 6 model solutions we 
selected 3 models as our 'favoured' ones. They have / = 1 
values for the 3 largest amplitude modes (J2, fi and /4). 
Fig. [3] shows the 6 models plotted by closed circles in 
the effective temperature - mass plane. The black square 
denotes the spectroscopic solution with its uncertainties. 
Our 3 favoured models are marked with open squares. 
As it can be seen in Fig. [3] two out of our 3 favoured 
models are too hot according to the 1 a limit in effective 
temperature. Therefore, the 0.645 Mq model seems to be 
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FIGURE 3. The 6 selected models in the T e ff - plane 
(closed circles). We denoted by black square the spectroscopic 
value with its uncertainties |6] |7]. Open squares mark our 
3 favoured models which give / = 1 values for the 3 largest 
amplitude modes. 



our best choice. 

We present further details of the frequency analyses, 
modeling and tests on amplitude variations of this star in 
Bognar, Paparo, Bradley & Bischoff-Kim J2]. 



INVESTIGATION OF LIGHT CURVE 
SEGMENTS 

We published our preliminary analysis of the first 
month's dataset in 2007 [3]. We investigated two seg- 
ments of the light curve which show alternating smaller 
and larger maxima. We found harmonic and lower- 
amplitude subharmonic peaks by the Fourier analyses of 
these segments. This feature of the spectra implies that 
the dynamics of the star is dominated by non-linear pro- 
cesses. 

Fontaine et al. |0] described KUV 02464+3239 as a 
photometric twin of the cool DAV star, GD 154. In the 
case of GD 154 Robinson et al. lfl4ll also detected har- 
monics and subharmonics in the power spectrum. The 
emergence of subharmonic peaks is similar to what we 
see in some non-linear dynamical systems which evolve 
toward chaos via a cascade of period doubling bifurca- 
tions. This result on the two segments made us to in- 
vestigate other parts of the light curve as well. However, 
we couldn't point out the presence of significant subhar- 
monic peaks but only of harmonics. This gave us the idea 
to investigate light curve segments from another point of 
view. 

We selected segments which show large amplitude 
variations. In these cases we can see easily the non- 
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FIGURE 4. Standard deviation values of the selected light 
curve segments. From top to bottom: segments with large am- 
plitude variations, segments showing beating of frequencies 
and small amplitude ones. Numbers of segments are inserted. 
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FIGURE 5. Normalized differential light curves of the 4 
largest amplitude segments. Numbers of segments are indicated 
in the upper-right corners of panels. 



sinusoidal shape of the cycles and we can find at least the 
second harmonic peak of the frequency which is domi- 
nant in the given segment. As control cases we selected 
light curve segments showing just small amplitude vari- 
ations and segments where we can determine two close 
frequencies causing beat in that part. 

We used the standard deviations of the light curve seg- 
ments to parametrize the groups. The standard deviation 
(SD) has two sources. The scatter of the measurement 
and the deviation from the mean light level caused by the 
pulsation. Since our observations are of high quality the 
dominant part of SD reflects the amplitude of the star's 
pulsation. Fig. [4]confirmes the correctness of our selec- 
tion criteria showing standard deviations of the selected 
light curve segments. The 3 groups are widely separated 
according to this parameter. The numbers in Fig. [4] are 
used to refer to a certain segment throughout the paper. 

Four segments with the largest amplitude light vari- 
ations are presented in Fig. [5] They have different rea- 
sons for being selected. Segment 5 analysed in Bognar et 
al. Ot] shows very explicit alternation of small and large 
amplitude cycles which denotes the presence of subhar- 
monic frequencies. In the case of Segment 1 1 we see an 
abrupt change in the peaks' amplitude from one cycle to 
another. Segments 40 and 41 seem to have a roughly si- 
nusoidal envelope suggesting strong beating of frequen- 
cies. In Segment 41 frequencies at ~87 and 104 c/d and 
their second harmonics are responsible for the envelope's 
shape. 

In the case of the large amplitude segments there are 
remarkable discrepancies between the observations and 
their fits. Neither the 6 accepted normal modes nor the 
extended list with 13 frequencies can give acceptable fits. 
SD seems to be a useful parameter again to show the 
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FIGURE 6. Investigation of the quality of fits in the case 
of large amplitude segments: difference values between the 
standard deviations of the original light curve parts and their 
fits. 



behaviour of segments using 6- and 13-frequency fits. In 
Fig. [6] we present the SD differences of the original light 
curve segments and their fits (solid and dashed lines: 6 
and 13-frequency solutions, respectively). Most of the 
segments show ~ 0.01 m SD improvements using more 
frequencies. The 13-frequency SD differences of most 
segments agree with the average SD difference of the 6- 
frequency fits for one of our control group, the beating 
segments. This level is ~ 0.014 m . 

Our examples presented in Fig. [5] have privileged po- 
sitions in the representation of SD differences, too. In 
Fig.|6]the 4 largest amplitude segments have the largest 
difference values. In the case of Segment 5 we see 
smaller improvement with a 13-frequency fit than for 
Segment 40 and 41. Considering the 13-frequency fit 
Segment 40 and 41 do not have privileged positions any- 



more. In some cases - at Segment 29, 26 and 37 - we get 
only slight improvements. The average of improvement 
using 13 frequencies instead of 6 is about 40%. 

CONCLUSIONS 

This short investigation of larger amplitude light curve 
segments suggests the importance of non-linear pro- 
cesses working in the star. Considering the shape of the 
light curve (with special regards to alternating smaller 
and larger amplitude cycles) and the FTs of segments 
showing large amplitude variations it is obvious that to 
understand the behaviour of this star we would need a 
non-linear treatment of its pulsation. 

Another result is that we found differences rather than 
similarities between the overall pulsational properties of 
GD 154 and KUV 02464+3239. The Fourier spectrum 
of GD 154 is dominated by only 3 modes and their har- 
monics and linear combinations lfl3ll . We found 6 modes 
in the case of KUV 02464+3239 and only one signifi- 
cant harmonic in the whole light curve's spectrum. We 
see clear presence of further harmonics only in segments 
of the light curve. The non-linear behaviour of KUV 
02464+3239 and similar white dwarf pulsators mean real 
challenge for theory and they are waiting for a solution. 
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